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We investigate electrical conduction through individual oligopeptide molecules thiol-bonded be- 
tween gold nanocontacts using ah initio and semi-empirical techniques. Our theory explains for the 
first time these molecules' experimentally observed current- volt age characteristics, including both 
the magnitude and rectification of the current, and uses no adjustable parameters. We identify the 
mechanism of the observed current rectification, and predict that it will result in negative differ- 
ential resistance at moderate biases. Our findings open the way to the realization of protein-based 
nanoelectronic devices. 



I. INTRODUCTION 

Living cells survive and function through their unique 
ability to manufacture a vast catalog of specific molecules 
repeatably and reliably, a capacity long recognized to 
hold promise for the design of nanoelectronic devices. 1 
Whereas much experimental and theoretical work has 
been devoted to understanding the electron transport 
properties of single nucleic acid molecules, relatively little 
has focused on the same qualities of other biomolecules. 
In particular, over the course of its life a cell produces 
far more polypeptide chains than DNA sequences, while 
amino acids, the building blocks of proteins, display five 
times more diversity than DNA bases. As such, the 
single-molecule electron transport characteristics of these 
protein fragments, and of entire proteins, is a field with 
great potential, which is only beginning to be explored. 

Recently, experimentalists have succeeded in address- 
ing the effects of primary protein structure, i.e. spe- 
cific amino acid sequences, on conductance^ These ex- 
periments, applying the STM-break junction (STMBJ) 
technique 4,5,6 to oligopeptide molecules, present an ex- 
citing new approach to biologically-based nanoelectronic 
devices: They offer a broad class of well defined experi- 
mental systems operating at length scales on which coher- 
ent quantum effects are dominant. Reference y, in par- 
ticular, demonstrated this promise with non-equilibrium 
results evincing a striking current rectification. 

The purpose of this article is to present the first theo- 
retical treatment of electron transport through these sys- 
tems. We explain the phenomena observed in these ex- 
periments, in particular the current rectification, by use 
of Landauer conductance theory^ and demonstrate good 
quantitative agreement with the experimental data. We 
show that the rectification is the result of resonance be- 
tween gold-molecule interfacial states partially localized 
on either side of the molecule (see Fig. [1]); these states 
are detuned from each other in equilibrium due to the 
oligopeptides' intrinsic asymmetry. Moreover, we predict 
that, at moderately higher voltages, a related mechanism 
will generate negative differential resistance (NDR), i.e., 
a decrease in current with increased bias voltage; this 
phenomenon has important device applications, includ- 
ing low-noise amplification, 9 high frequency oscillators^ 
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FIG. 1: (Color online) Anticrossing of generalized eigenstates 
R and L of the simplified extended molecule discussed in the 
text for cysteamine-Gly-Cys bonded to two hollow sites in the 
Au(lll) leads. In the full system, the anticrossing's abrupt- 
ness is tempered by broadening from the leads. The insets 
show the states at specific bias voltages. 8 



analog-to-digital conversion^ and digital logic. 12 Thus, 
our findings open the way to the realization of practi- 
cal oligopeptide, and ultimately protein-based, nanoelec- 
tronic devices. 

We organize the paper as follows. Section [III presents 
the theoretical formalism we use to model charge trans- 
port through metal-oligopeptide-metal junctions. In 
Section we discuss our results for these systems, in- 
cluding determination of the experimental geometry, ex- 
planation of their rectifying properties, and prediction of 
NDR. Section HVl summarizes our main conclusions. 
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II. MODEL 

To determine the molecular geometry in the presence 
of the gold leads, we apply a series of density functional 
theory ab initio calculations to relax the molecule be- 
tween two fixed gold atoms At each step, the dis- 
tance between the two atoms is increased slightly, and 
the geometry relaxed again, to simulate the stretching the 
molecule undergoes in an STMBJ experiment. For calcu- 
lations reported here, the molecules have been stretched 
to their greatest length before breakage occurs, although 
we find that less stretched molecules have essentially the 
same conductance features and order of magnitude. In 
the case of hollow binding sites, we replace the gold atom 
with the three nearest-neighbor gold atoms, and relax the 
molecule one final time. 

To determine conductance and current, we first parti- 
tion the metal-oligopeptide-metal system into left and 
right macroscopic gold electrodes, and an extended 
molecule containing the central organic molecule as well 
as 50-100 atoms of each Au(lll) contact. 14 The Hamilto- 
nian of the full system is therefore the sum of two terms: 
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H mo i describes the isolated extended molecule, and 
Hieads does the same for the macroscopic electrodes, as 
well as containing the quantum mechanical coupling be- 
tween them and the extended molecule. 

We employ a semi-empirical model for H mo i; such 
models are known to be highly effective in predict- 
ing, reproducing, and explaining experimental results for 
metal-organic molecule-metal junctions (e.g., Ref. [TBI, 
Ref. M and the results of Refs. O andi, Ref. [H). The 
Hamiltonian of the full extended molecule is given by 
extended Hiickel theory^^ 2 ^ 2 ^ 3 .^ 

Hieads provides both decoherent broadening and elastic 
scattering to the discrete states of the extended molecule. 
To model this, we assign semi-infinite, one-dimensional 
leads 25 to each gold orbital not immediately adjacent 
to the oligopeptide; their number is equal to the coor- 
dination number of gold minus the number of nearest- 
neighbor atoms included in the extended molecule. We 
take the intersite hopping within each lead, as well as 
the hopping between lead and orbital, to be 5eV. These 
choices are made to quantitatively reproduce the quan- 
tized conductanc e) 26 ! 27 of a gold quantum point contact, 
and are not fitted to any molecular electronics experi- 
ment. 

The open system described by H requires an infinite 
Fock space. To solve the full molecule-plus- leads system, 
then, we make use of Dyson's Equation, which gives the 
retarded Green function 
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Here 
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is the Green function for the extended molecule onl yr^ 
and the retarded self-energy T,(E) is given by Ref. [25|. 
S is the atomic orbital overlap matrix of the extended 
molecule. Equation ([2]) thus represents an exact method 
to incorporate the effects of Hieads into transport. 

We model electrostatic effects, including the potential 
imposed by the contacts as well as screening by the free 
ionic charges in the electrolyte solution surrounding the 
molecule, with a linearized Debye-Hiickel theory* 2 ^ In 
such a theory, Poisson's equation for the total electric 
potential <E>(r) becomes 
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where n = .1M is the concentration of the solution, — e 
is the charge of an electron, e is the dielectric constant 
of water, ks is the Boltzmann constant, and T is room 
temperature^ 

We solve Eq. (j4j) for parallel-plate boundary condi- 
tions, where the plates are taken to lie at the metal- 
molecule interfaces, and take $(r) constant within the 
metal nanoclusters themselves. In the electrostatic gauge 
for which the voltages of the left and right lead add to 
zero, the solution of Eq. (jlj) in the region containing the 
molecule is 
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where z = is one boundary surface, z = d the other, V 
th e applied b ias, and the Debye screening length is A = 
Y / /c J eTe/2ne 2 . For the experimental systems of Refs. 2 
and [3|, A = 0.97nm; since this length is greater than 
d/2 for all the metal-molecule-metal junctions studied 
here, the effect of screening in these systems is slight, 
and our results are robust against reasonable changes in 
the experimentally determined parameters T, e, <i, and n. 
It should also be noted that, although Eq. ([5]) is written 
as an example in a particular gauge, the results of our 
calculations are gauge invariant. 

We use a simple ansatz to include the electrostatic ef- 
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where, is the position of the i th atomic orbital. Thus, 
each diagonal element is gated in the usual way, and 
off-diagonal elements are modified according to the non- 
orthogonality of the basis states, as is typical in a Hiickel 
theory. In addition, of course, the site energies of each 
lead and metallic nanocluster shift uniformly with their 
respective metallic contact's electrostatic potential. All 
the quantities in Eq. ([2]) thus become functions of the 
applied bias V, without the use of free parameters. 

The Landauer formalism 7 - yields the conductance and 
current through the molecule. The probability for an 
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TABLE I: Comparison of the conductances found in our model 8 (in units of go) for various bonding geometries with those of 
the STMBJ experiments.—^ In the column headings, the first entry indicates the bonding geometry on the amine side of the 
molecule, the second the other. Only hollow bonding at both interfaces is compatible with experiment. 



molecule 


experimental 


hollow- hollow 


hollow-top 


top-hollow 


top-top 


cysteamine-Cys 


1.8xl0" 4 


5.4xl0" 5 


1.5xl0" 6 


9.6xl0" 7 


5.8xl0" 9 


cy st eamine- Gly- Cy s 


4.2xl0" 6 


4.3xl0" 6 


5.3xl(T 8 


3.7xl0" 8 


1.5xl0" 9 


Cys-Gly-Cys 


5.3xlCT 6 


4.2xl0" 6 


9.3xl0" 8 


l.OxlO" 7 


1.7xl0" 9 


cyst eamine- Gly- Gly- Cys 


5.0xl0" 7 


4.1xl0" 8 


7.3xl0" 10 


1.9xl0" 10 


2.2xl0" n 



electron of energy E to move from the left set of leads to 
the right is 



a e left 
j3£ right 
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where 



r (a) (E,V) = -21m \^ a) (E,V) 
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is the broadening function for the term of E from lead 
a. The zero-bias conductance of the metal-molecule- 
metal junction is g = go^(^F,0), with the conductance 
quantum go = 2e 2 /h and the bulk gold Fermi energy Ep 
evaluated for a 459-atom spherical gold nanocluster. The 
current is 7 



I(V) 



dET(E, V) [f(E-n R )-f(E-VL) 
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where f(E — Hr,l) is the Fermi function of the leads 
on the right (left) side of the extended molecule, and 
Hr,l{V) is the corresponding electrochemical potential. 



III. RESULTS AND DISCUSSION 
A. Conductance and Bonding Geometry 

We studied the conductance and current of each 
oligopeptide molecule measured experimentally in Refs.lH 
and [3, varying several parameters of the experimental ge- 
ometry in each case. These included the sulfur-to-sulfur 
distance as the molecule was stretched, the presence of 
additional molecules of the self-assembled monolayers, 
and bond geometry; only the last was found to have an 
impact on the observables of more than a factor of order 
unity. Figure [2] demonstrates this dependence, consid- 
ering all four possibilities involving the two most likely 
bonding geometries of the sulfur^ either directly to a 
single gold atom on top of the Au(lll) substrate, or 
above a triangular hollow. Table U compares the con- 
ductances obtained from our modei^ with experiment. 

The experimental data are consistent only with the 
model calculations for a hollow bonding geometry at 
both ends of each oligopeptide molecule. Since this 
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FIG. 2: Equilibrium transmission functions 8 T(E, 0) for var- 
ious bonding geometries for each of the four oligopeptide 
molecules studied experimentally in Refs.Q and0. Solid lines 
indicate hollow-site bonding geometry on both sides. Dot- 
ted: top-site at both metal- molecule interfaces. Long dashes: 
hollow-site bonding for the sulfur on the amine-terminated 
side of the molecule and top-site bonding for the other. Short 
dashes: vice versa. Experimental conductances agree well 
with the solid curves in all cases. For the hollow-hollow ge- 
ometry, the energies of the two molecular orbitals responsible 
for rectification are indicated with crosses. Inset: Structural 
formulae of the gas-phase molecules. When bonded to gold, 
the terminating hydrogen atoms are lost. 



model is known to produce quantitatively correct con- 
ductance in other gold-thiol systems ^ 15 i 16 i 18 i 33 we con- 
clude that this is indeed the experimental geometry, and 
that the orders-of-magnitude smaller conductances as- 
sociated with top-bound geometries and with molecules 
bound via monatomic gold chains^ form part of the low- 



conductance noise in the experimental histograms^ In 
this way, the peaks of the conductance histograms of the 
STMBJ technique correspond to a specific geometry for 
the metal-molecule-metal junction, despite sulfur's abil- 
ity to bond to a gold substrate in several geometries. 

Further evidence that the single oligopeptide molecules 
of the experiments are bound to hollow substrate sites 
at each end is found in the behavior of conductance as 
a function of molecular stretching. Using DFT relax- 
ation at progressively larger lead- lead distances, we find 
that internal degrees of freedom within the oligopeptide 
molecule itself accommodate only about one third of the 
total stretching induced by pulling the leads apart. The 
remaining stretching is allowed only by repositioning of 
the entire molecule between the two metal contacts, i.e. 
by changes of the Au-S-C bond angles and the dihedral 
angles involving both gold contacts and the molecule. 
This finding is consistent with the persistence lengths of 
polypeptide chains, typically several nanometers or more. 

Indeed, the tetrahedral sp 3 bonding of the hollow- 
bound geometry is much less flexible, and its symmetry 
more restrictive, than the linear sp bond of the top-bound 
geometry. Therefore, when the hollow-hollow geometry 
is stretched, the DFT geometry optimizations only allow 
stretching of .2A or less (which affects T(E) by 10% or 
less) before the molecular wire detaches from one elec- 
trode. By contrast, we find stretching distances as large 
as lA in top-bound geometries, and the conductance of 
unstretched top-bound geometries can be as much as 3-4 
times larger than that of the fully stretched configura- 
tions reported in Table H as shown in Fig. OH Since the 
experimental data^ 3 - exhibit no strong stretching depen- 
dence of conductance within a single plateau, we conclude 
that the experimental results correspond to the hollow- 
hollow geometry^ in agreement with the conclusion 
reached above based on the magnitude of conductance^ 

Kriiger et al. have suggested^ that when a molecu- 
lar wire thiol-bonded to gold contacts is stretched and 
eventually breaks, the cleavage does not occur at the 
gold-thiol interface, but instead that a gold atom, re- 
maining attached to the thiol group, is pulled from the 
contact. During this process, the system passes through 
several intermediate, temporary configurations, including 
structures resembling the top-bound geometry discussed 
above. While a complete analysis of this geometric evo- 
lution and accompanying conductance changes is beyond 
the scope of the present work, we note that the results 
shown in Table [J indicate that the transition from hollow- 
hollow to a top-bound geometry at either end is associ- 
ated with a drop in conductance of more than an order 
of magnitude in each oligopeptide molecule. Thus, in 
a stretching experiment, when the top-bound geometry 
forms, the conductance plateau ends and the molecule 
appears to have simply separated from one or both con- 
tacts. 

Furthermore, during stretching, temporary intermedi- 
ate configurations such as those in Ref. [34] may form, 
which in general are less conducting than the original 
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FIG. 3: (Color online) (a) Equilibrium transmission functions 
T(E, 0) for cysteamine-Gly-Cys in the top-top configuration, 
for different distances between the two gold leads. Conduc- 
tances of the metal-molecule-metal junction, shown in inset 
as a function of tip-tip distance, vary approximately linearly 
over a range of about lA, beyond which no stable molecular 
structure bridges the two electrodes. Thus, the amount of 
stretching can affect conductance by a measurable amount, 
though not sufficiently to reconcile single- Au bonding geome- 
tries with the experimental results. For the hollow-hollow 
configuration, on the other hand, this range is much smaller, 
and so almost no variation in conductance is allowed, (b) and 

(c) show the L and R states, respectively, for the top-top sim- 
plified extended molecule with a tip-to-tip distance of 14. 2A. 

(d) and (e) show the same for a tip-to-tip distance of 15.2 A. 
Most of the stretching is accommodated by the bond angles 
at the gold-molecule interfaces, so that the principal change 
in (a) is peak widths. 



configuration. One probable intermediate bonding ge- 
ometry has the sulfur directly above a single gold atom 
embedded within the top atomic layer of the contact; 
when one end of the molecule is bound in such a way 
to the gold, and the other is in the hollow configuration, 
we find the molecular junction has a conductance of 1/4 
that of the hollow-hollow configuration. Indeed, weak 
features of lower conductance are seen in some of the ex- 
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perimental conductance vs stretching curves of Ref. 0. 
They are typically much shorter than the main plateau, 
and occur at its end, consistent with configurations that 
are only meta-stable and occasionally populated via the 
stretching mechanism. 

To understand the transmission curves, it is helpful 
to think in terms of a smaller extended molecule which 
includes only the two to six gold atoms to which the sul- 
furs bond. Transmission states near Ep correspond to 
gold-sulfur interfacial orbitals of this simplified extended 
molecule. Because there are two metal-molecule inter- 
faces, these states typically occur in pairs with similar en- 
ergies. When the two states are detuned from each other, 
they are largely localized at either end of the molecule; 
when in resonance, they mix into bonding and antibond- 
ing combinations (see Fig.[T]). For a symmetric extended 
molecule at zero bias, the pairs are always in resonance; 
oligopeptides, however, are intrinsically asymmetric, so 
interfacial states can be detuned from each other and 
localized at either end, even in equilibrium. 

Prominent peaks in the transmission spectra (Fig. [2j) 
occur near the energies of molecular orbitals which bridge 
the left and right leads. Thus, the resonances, and conse- 
quent hybridization, of levels otherwise localized on op- 
posite ends of the extended molecule play a dominant role 
in determining electron transport. For example, in all 
molecules and bonding configurations studied here, the 
two levels immediately below Ep are paired interfacial 
states, which we call R (localized at the gold-molecule 
interface near the amine terminus of the molecule) and 
L (localized at the other interface). The narrow trans- 
mission peaks associated with their bonding and anti- 
bonding combinations are visible leV below Ep in the 
double on-top (dotted) curves of Fig. [2j Because the 
overlap of the gold 6s and sulfur 3p orbital is strongly 
angle dependant, 37 these peaks each become broadened 
when their respective bonding geometries are changed 
to hollow. Furthermore, the hollow-bonded versions of 
R are more sensitive to the asymmetric features of the 
oligopeptide molecules, resulting in a significant detuning 
from L (shown by the crosses in Fig. [2]). The transmission 
peak associated with these two states is thus also severely 
broadened, becoming a smooth shoulder in T(E,0). 



B. Current Rectification 

The movement and change under bias of R and L, 
depicted in Fig. [TJ give rise to the rectification observed 
by Xiao et al in Ref. 13. Positive bias raises (lowers) 
the electrochemical potential of atomic orbitals on the 
amine (non-amine) side of the molecule; consequently, 
under negative bias the energies of R and L move closer 
together. Near resonance, R and L hybridize into whole- 
molecule states, allowing significant charge transport. 

Figured! shows the movement under bias of several gen- 
eralized eigenenergies of the simplified extended molecule 
for hollow-hollow bonded cysteamine-Gly-Cys. 8 All of 
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FIG. 4: Movement with applied bias of the general- 
ized eigenvalues^ of the six-gold-atom simplified extended 
molecule for cysteamine-Gly-Cys in the hollow-hollow config- 
uration. Some lines are doubled by degrees of freedom in the 
gold contacts; when more gold atoms are included, the molec- 
ular orbitals split further, eventually developing into the con- 
tinuous density of states of an infinite system. Note that the 
apparent intersections are actually avoided crossings. In equi- 
librium, states R and L are detuned by .4eV due to molecular 
asymmetry. They come into resonance only under negative 
bias, causing current rectification. Without these two states, 
the location of anticrossings, and thus the I-V characteristic, 
would be largely symmetric. The dashed lines represent the 
electrochemical potentials of the right and left electrodes, /hr 
and NDR is the results of increasing bias beyond the last 
anticrossing in a group; an example is circled. 



the states in the region around Ep occur in right/left 
pairs like R and L, but only R and L are detuned from 
each other in equilibrium. Without R and L, the loca- 
tions of anticrossings are nearly symmetric about V = 0. 
Rectification is a result of the splitting between these two 
states, itself a direct consequence of oligopeptide asym- 
metry. 

Figure [5] gives the calculated current- volt age character- 
istic for each molecule^ showing this rectification. While 
the resonance depicted in Fig. [1] is comparatively narrow 
and outside the energy window between \ip and recti- 
fication extends beyond 2V in the full systems due to the 
strong hollow-site broadening of R and L. For the three 
molecules for which Xiao et al. made I-V measurements 
(shown in inset), the calculated rectification is in excel- 
lent qualitative agreement with experimental data. We 
find a current ratio of about 150% for cysteamine-Gly- 
Cys and Cys-Gly-Cys, similar to Ref.0. Cysteamine-Gly- 
Gly-Cys exhibits less rectification, both in our model and 
the experimental data: This is due to its much smaller 
detuning between R and L at zero bias. 



C. Negative Differential Resistance 

Negative differential resistance (dl/dV < 0) can result 
from applied bias destroying a resonance between trans- 
mitting interfacial states j 38 i 39 As the bias increases past 
an anticrossing in the metal-oligopeptide-metal systems, 
the two interfacial states are taken further apart in en- 
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FIG. 5: Calculated current- volt age characteristics of four 
single-molecule oligopeptide devices.— The scale and asym- 
metry of the current are comparable to those found in the 
experimental studies of Ref. 3, reproduced here in inset. 



ergy. They therefore mix less well, and become localized 
at the individual metal-molecule interfaces. As the bias 
voltage increases, then, electron transport through these 
states "shuts off" , creating the potential for NDR. 

Such effects are masked while the current is still small: 
The dominant effect is the increasing current from a se- 
ries of crossings and states entering the window of inte- 
gration. Nevertheless, once the voltage is large enough 
that an entire neighborhood of interfacial resonances has 
been encompassed, further applied bias yields NDR. One 
such crossing is indicated by the circle in Fig. [4j its I- 
V characteristic and relevant transmission properties are 
shown in Fig. [6l As V becomes more negative, the peak 
in T(E) associated with the nearly degenerate hybridized 
states diminishes. After the final current step, due to a 
mid- molecule orbital (★ in Fig. 2j) entering the window of 
integration, NDR results. The bias voltages required to 
observe this important phenomenon are only moderately 
larger than those already achieved experimentally^ 

The magnitude of NDR is determined by the partici- 
pating levels' electrostatic and quantum mechanical cou- 
plings to the leads. While the first is necessarily high 
for interfacial levels (hence the saturated slope of most 
lines in Fig. the quantum mechanical coupling, or 
broadening, of an interfacial level depends mainly on the 
bond geometry at the metal-molecule interfaces. For hol- 
low bonding, for example, interfacial states are very well 
broadened, as can be seen from the hollow-hollow (solid) 
curves of Fig. [2j On the other hand, if one or both of 
the interfaces is in the top-bound geometry, the related 
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FIG. 6: I-V characteristic for the hollow- hollow gold- 
cyst eamine-Gly-Cys-gold system^ in the medium-bias regime, 
showing NDR. Other oligopeptide molecules have qualita- 
tively similar characteristics, (inset) Corresponding T(E, V) 
curves. NDR corresponds to a gradual decrease in the trans- 
mission peak's strength as voltage becomes more negative, 
due to departure from a resonant crossing of interfacial states 
(the circle in Fig. |4}. 



NDR is stronger. 

Like their rectifying properties, the NDR that we 
predict in oligopeptides is a direct result of the domi- 
nance of resonant crossings of interfacial states in charge 
transport. Since this is a property of all thiol- (i.e. 
cysteamine- or cysteine-) bonded oligopeptide molecules, 
both effects are robust, and observable in the entire class 
of such molecules. Because of the limitless customizabil- 
ity and simple fabrication of oligopeptides, the effects can 
be tailored to meet a wide variety of device needs. 



IV. CONCLUSIONS 

In summary, we have developed a theoretical model 
which explains all of the experimental data on charge 
transport through oligopeptide molecules. We obtain 
quantitative agreement with experiment by using a com- 
bination of appropriate, standard theoretical models 
(DFT for the molecular geometry; the linearized Debye- 
Hiickel model for the potential profile, including the 
electrostatic screening effect of the electrolyte solution; 
and the extended Hiickel model of Hoffman and co- 
worker o 19 i 20 for the electronic structure) and careful con- 
sideration of the experimental systems. We use no fitting 
of any kind to any experimental data involving electrical 
conduction through molecules. Quantitative agreement 
without the use of fitting parameters allows us to extract 
detailed chemical and physical information from the ex- 
periments, such as the conclusion that both ends of the 
molecule are bound to hollow sites on the Au(lll) sur- 
faces. 

We understand the observed current rectification of 
oligopeptides as due to a hybridization of interfacial 
states lying at either end of the molecule, naturally de- 
tuned from each other in equilibrium due to the intrin- 
sic asymmetry of peptide chains. When bias is applied 
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in one direction, these states are brought into resonance 
with each other, forming two whole- molecule states which 
allow a stronger current to flow; such a resonance does 
not occur when the molecule is biased in the opposite 
direction. Our calculated current rectifying ratios agree 
well with those measured in experiment, again without 
the use of fitting parameters. 

Finally, our model predicts NDR will appear at moder- 
ately higher bias voltages than have been explored exper- 
imentally in Refs. and 0. This phenomenon is under- 
stood as the natural counterpart to the rectifying physics: 
As bias is increased past the anticrossing of the two in- 
terfacial states, the states once again become localized at 
the metal- molecule interfaces, and current flow through 
these channels decreases. Since NDR has many atten- 
dant device applications, and oligopeptides offer a unique 
combination of customizability and scalability, it is our 
hope that this work will inspire much further study of 



their electron transport properties. 
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